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With a sufficiently high number of fundamental fermionic flavours present, Yang-Mills theory 
develops an infrared fixed point and becomes (quasi-)conformal in nature. The range of flavour 
numbers for which this occurs defines the conformal window, the lower limit of which has yet to 
be determined. We studied the phase diagram of SU(3) Yang-Mills theory with twelve flavours 
of staggered fermions. Here we present evidence for the occurrence of a bulk transition towards 
a chirally broken phase and the existence of a Coulomb phase on the symmetric side of this 
transition, using results from the measurements of the chiral condensate and spectrum, leading to 
the determination of a positive sign of the beta function. Assuming the validity of the Appelquist- 
Miransky-Yamawaki scenario, this implies the existence of a conformal window that comprises 
the theory under investigation. 
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There has been a recent surge in interest in the physics of gauge theories with a relatively large 
number of fermionic degrees of freedom. Research is ongoing on theories with different gauge 
groups and matter content in several representations [1,2]. The common factor of interest for all of 
these is the possible appearance of so called Banks-Zaks fixed points, first described in two seminal 
papers of the late seventies and early eighties [3, 4]. The central observation of these papers was the 
appearance of a second zero before the loss of asymptotic freedom in the two-loop beta function of 
an SU(3) gauge theory with Nf massless fermions in the fundamental representation. 

Restricting ourselves to the latter case, the beta function to second order in perturbation theory, 
the highest order to which all coefficients are universal, is given by 



This equation has a root for Nf > 8.05 at g* 2 = —bo/bi, causing, at this perturbative level, the 
appearance of an infrared fixed point up until Nj = 16.5 when asymptotice freedom is lost. Physi- 
cally, a IR fixed point implies a phase in which chiral symmetry is never spontaneously broken, as 
the coupling reaches an asymptotic value at infra-red scales and the theory turns quasi-conformal. 

The actual value of the critical number of flavours Nj is modified by the non-perturbative 
dynamics of chiral symmetry breaking and its restoration [5], necessitating an approach beyond 
perturbation theory. Analytical approximations have been done, generally favoring a number larger 
than the perturbative prediction [6-8], but with no consensus on the actual value. Lattice studies 
have concluded that Nf = 8 lies within the hadronic phase of QCD [9- 13], but the case of Nf = 
12 is debated. A study based on the Schrodinger functional [9, 10] has concluded that Nf = 12 
should already be in the conformal window, but other groups have reported indications of a regular 
hadronic phase and spontaneous chiral symmetry breaking instead [12-14]. We here report on an 
investigation of the nature of SU(3) Yang-Mills theory with twelve flavours of quarks, employing 
the physics of phase transitions, to appear in print shortly [15]. 

1. Setup and implementation 

In the setup of this study, we took guidance from the scenario painted in [5], a diagramatic 
depiction of which is given in figure 1. For a given value of the flavour number Nf > Nj and in 
the limit of massless quarks, decreasing j8 will cause one to cross the critical value and go from 
the chirally symmetric (S) and asymptotically free phase to a phase that is still symmetric, but not 
asymptotically free. The green dashed line connects these roots of the beta function. The numerical 
values for which these occur are action dependent, but since the presence of such a fixed point is 
physically observable, it will appear for some value at any choice of the action. Note, however, that 
the beta function is not universal away from fixed points with diverging correlation lengths and one 
can therefore not exclude the appearance of additional, spurious fixed points [16]. At even smaller 
values of /3 , the lattice regularized theory could in general exhibit a transition to a strongly coupled 
chirally asymmetric (A) phase, driven purely by the value of coupling constant itself. 
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Figure 1: Diagram of the phases of an 51/(3) gauge theory encountered through varying the flavor content 
Nf in the fundamental representation and the strength of the coupling g. Theories for Nf < Ni are QCD-like, 
while for iVS > Nf > develop a conformal phase. S and A refer to chirally symmetric and asymmetric, 
respectively. The dot-dashed (red) line is the line of lattice bulk transitions. We observed a bulk transition 
at Nf = 12 and also 16 (dots). The dashed (green) line is the conformal critical line, qualitatively indicating 
the location of the IRFP's. The line at Nf = NJ indicates a conformal phase transition with abrupt change 
of the spectrum, predicted in reference [5]. The non-perturbative beta function on the conformal side is also 
shown. 



Outside the conformal window, in contrast, the theory will break chiral symmetry at zero 
temperature, but exhibit a thermal phase transition in the continuum to a high temperature chirally 
symmetric phase. Thus, as was argued in [11], observing such a transition implies the absence of 
an infra-red fixed point. To discriminate between a true thermal transition of the hadronic phase 
and the aformentioned bulk transition of the quasi-conformal phase, we observed the scaling of the 
transition with the physical temperature, given by T = (a(f5)N T ) i . The two types of transition 
are distinguished by their zero temperature, or infinite N x limit. If we obtain lim^^oo jSr — > °°, we 
are in the hadronic phase, while lim/v^co j3r = fi c implies a bulk transition for a finite value of J3 C . 

The bulk transition is a lattice artifact, one that specific forms of improvement may push 
down to arbitrarily low values of the lattice coupling and possibly have vanish altogether. Its 
potential appearance is relevant in that it would exclude its thermal counterpart. However, it is not 
inconceivable that one observes a transition between two different chirally broken phases, leaving 
room for a weak thermal transition at some smaller value of the coupling. This scenario can be 
excluded primarily by checking directly for chiral symmetry restoration. In addition, the presence 
of a Coulomb phase, signalled by a positive value of the beta function, is a necessary and sufficient 
condition for the presence of a Banks-Zaks fixed point in the discretized system. 

Simulations were performed with twelve dynamical flavours of fermions in the form of three 
unrooted staggered fermions. A tree level Symanzik improved gauge action was used to suppress 
lattice artifacts. Fermions were Naik improved for an overall 6 {a) improvement that should be 
valid in any phase of the theory. Allowing for the appearance of phases outside of regular hadronic 
QCD, other types of improvement may introduce unphysical terms in the continuum lattice ac- 
tion. Specifically, the Lepage-Mackenzie tadpole improvement scheme, an obvious choice when 
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Figure 2: The bulk transition in the chiral condensate for am q = 0.05 at several values of the temporal 
lattice extent N T . The location of the transition is identical, in spite of the physical temperatures differing 
by substantially. Simulation errors are within symbol sizes. 



attempting to observe thermal scaling for a phase transition, but derived from a perturbative expan- 
sion [17], may introduce unphysical effects when used in the absence of asymptotic freedom. Taste 
breaking effects were estimated by comparing different pion sources and found to be of the order 
of 40%. By increasing pion masses, the effect of taste breaking contributions should be to partially 
freeze chiral degrees of freedom and in this way effectively lower the number of fermionic degrees 
of freedom. Any conclusions drawn from the measurements presented here will therefore specifi- 
cally apply to twelve staggered degrees of freedom at finite mass, with a non-trivial continuum and 
chiral limit. However, the shielding effect of the fermions responsible for any fixed point can only 
increase in those limits and if indications of its existence are indeed found, they will carry through. 

Runs were done on several volumes, with explicit checks on the size of finite volume effects 
for the high precision measurements needed to perform the chiral extrapolation by comparing the 
results of three volumes and finding agreement within statistical errors for the largest two. Final 
fits were based on 24 4 volumes. Spectral measurements were done on 16 3 x 24 volumes, intro- 
ducing a moderate amount of finite volume corrections to the values obtained that tend to increase 
and equalize the obtained masses. Those effects, as long as they are not dominating, will wash 
out the signal obtained for the sign of the beta function, but should not be expected to change it 
qualitatively. 



2. Results 



The size of the chiral condensate (yy) was determined over a series of values of the lattice 
coupling constant j3 (see fig. 2) and found to be a smoothly varying quantity for all measurements. 
In the range between /3 = 3.0 and j8 = 3.5, the condensate experienced a rapid acceleration in its 
decline, that tapered off again towards the region of weaker coupling. As observed above, such 
a transition can occur for actions within and outside of the conformal window, but is driven by a 
different quantity for either case. To see wether the physical temperature was the relevant factor, 
measurements were repeated for values of the lattice extent in the temporal direction N r differing 
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Figure 3: Extrapolation to the chiral limit of the chiral condensate for j3 = 3.9 with a linear formula 
(reduced % 2 = 3.81), giving a small amount of residual chiral symmetry breaking, and a power law (reduced 
X 2 = 0.86), demonstrating a strong preference for the chirally symmetric scenario with a non-zero anomalous 
dimension. 

by more than a factor of two. While the measurements on the weak coupling side of the transition 
show clear thermal sensitivity, we observed no associated shift in the position of the transition 
itself. This implies that lim^ T ^oo fir remains finite and the transition is of bulk nature. 

To check for the restoration of chiral symmetry, an extrapolation to the chiral limit was per- 
formed. The dependence was approximately linear, but a slight systematic deviation of this be- 
haviour could be observed. Most notably, higher masses for coupling showed clear upward de- 
viations from linear behaviour close to the transition. This behaviour was reported earlier by 
Damgaard et al. [16], and is explained by both a shift in the position of the transition to weaker cou- 
pling and a broadening of the transition region as masses become heavier, implying one observes 
remnants of the transition itself. Away from the transition these deviations rapidly weakened and 
the data developed a convex deviation from linearity. As these measurements were performed in 
an interacting system, one should generally allow for a non-zero anomalous dimension, as was, 
for example, shown for the chiral transition in the case of strongly coupling QED [18]. Additional 
high precision measurements were performed after the meeting, to be reported upon elsewhere, 
to perform a quantitatively reliable chiral extrapolation. Figure 3 displays an extrapolation to the 
chiral limit for a single value of the coupling of j3 = 3.9 using, both with a linear form and allowing 
for an anomalous dimension. The difference between both forms is small for the available masses, 
but the data prefer a scenario without spontaneous chiral symmetry breaking with a fitted exponent 
of 0.965(1). 

Having observed a bulk transition to a phase that is compatible with the restoration of chiral 
symmetry, we tried to qualify the weakly coupling phase further. One can attempt to distinguish 
a very weakly coupling hadronic phase from a Coulomb phase directly, by determining the sign 
of the beta function. The beta function itself is a scheme dependent quantity, even when it is de- 
fined non-perturbatively, but the physical nature of the fixed point should make the sign of the beta 
function a rather robust feature. We set up such a scheme by following reference [16], measur- 
ing the spectrum of light mesonic excitation at different values of the coupling constant and bare 
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Figure 4: The measured ratio of m vt and m ps as a function of the bare quark mass, for a series of values of 
the coupling constant /3 . Identical values of this ratio define theories extrapolating to the continuum limit, 
sampled at different scales. The trend is towards higher values of j3 with higher masses matching with lower 
values of beta and lower masses, indicating a decrease of the interaction strength for coarser lattices. 



quark mass. Those coordinates in the two dimensional parameter space connected by renormal- 
ization group transformations should produce identical physical values, modulo a scaling factor. 
Subsequently, any dimensionless ratio of physical observables should be a constant, assuming that 
the system can be described by one parameter beta functions. Choosing the ratio of the ground 
state excitations in the pseudoscalar and vector channels arbitrarily and assuming a locally smooth 
behaviour of the function, the argument can be inverted to determine sets of parameters defining 
identical physics at different length scales. 

The spread in values for the ratio is smaller than that observed for the case of Nf = 16 in [16], 
but the overall trend is identical in that theories at weaker coupling and lower mass match to those 
at stronger coupling and higher mass. Because the signal-to-noise ratio is lowered, the result is 
not as outspoken. One way of improving this is by interpolating the values of the ratio, effectively 
combining statistics for the measurements. Such an extended analysis will be published [15]. 
Comparing the ratios in 4, the phase after the bulk transition is identified as a Coulomb phase. 
Assuming lattice artifacts are not dominating and under the assumption of continuity (which can 
be checked by scanning the parameter space towards zero coupling), this phase should connect 
smoothly to the perturbative function as the interaction weakens over larger distances. Since its 
sign is known to be negative perturbatively an intermediate zero has to be present. 
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